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CLOCK SKEW MEASURING APPARATUS AND METHOD 



BACKGROUND OF THE INVENTION 



5 1. Field of the Invention 

The present invention relates to a clock skew measuring 
apparatus and a clock skew measuring method for measuring skews 
between a plurality of on-chip clock signals obtained by- 
distribution by a clock distribution circuit on a semiconductor 
10 integrated circuit chip. 

2. Description of the Related Art 

According to conventional techniques, an on-chip clock skew 
is statistically estimated by bringing out clock signals to be 

15 measured into the outside of a chip at the same time, as shown 
in Fig. 1, using a time interval analyzer or frequency counter. 
The time interval analyzer measures a timing difference of a 
zero-crossing point between the clock signal to be measured and 
a reference clock signal so as to measure fluctuation thereof in 

20 accordance with histogram analysis. An exemplary clock skew 
measuring technique using the time interval analyzer is disclosed 
in Wavecrest Corp., Jitter Analysis Clock Solutions, 1998, for 
example . 

25 The conventional clock skew measuring method, however, 

requires a plurality of high-frequency clock output pins that may 
increase a cost of the apparatus, in order to simultaneously bring 
out a plurality of clock signals to be measured to the outside 
of the chip, thereby increasing the cost of the measurement. 

30 Moreover, the number of pins of the chip is limited. Thus, in 
a large-scale semiconductor integrated circuit, the skew 
measurement is possible only for a very small number of the 
distributed clocks, so that the clock skew for the entire circuit 
cannot be obtained precisely. Therefore, a new clock skew 
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measuring method is required in order to enable the high-precision 
control of the skews between the on-chip clock signals. 

It is an object of the present invention to provide a clock 
5 skew measuring apparatus and a clock skew measuring method which 
can efficiently estimate the skews between the on-chip clock 
signals . 

SUMMARY OF THE INVENTION 

10 

Therefore, it is an object of the present invention to provide 
clock skew measuring apparatus and method which overcomes the above 
issues in the related art . This obj ect is achieved by combinations 
described in the independent claims. The dependent claims define 
15 further advantageous and exemplary combinations of the present 
invention. 

According to the first aspect of the present invention, a 
clock skew measuring apparatus for measuring a clock skew between 

20 a plurality of clock signals to be measured in a device under test, 
comprises: a clock signal selecting element operable to receive 
the plurality of clock signals to be measured and to output the 
plurality of clock signals to be measured by selecting one of the 
plurality of clock signals to be measured one by one; and a clock 

25 skew estimator operable to receive a reference signal input to 
the device under test and the plurality of clock signals to be 
measured output from the clock signal selecting element and to 
obtain the clock skew estimates between the plurality of clock 
signals to be measured by measuring a timing difference between 

30 the received reference signal and each of the plurality of clock 
signals to be measured. 

The clock skew measuring apparatus may further comprise: 
a plurality of buffers operable to respectively supply the 
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plurality of clock signals to be measured to the clock signal 
selecting element; and a controller operable to control whether 
or not each of the plurality of buffers supplies a corresponding 
one of the plurality of clock signals to be measured to the clock 
5 signal selecting element. 

The clock skew estimator may measure a deterministic 
component of the clock skew between the plurality of clock signals 
to be measured. 

10 

The clock skew estimator may measure a random component of 
the clock skew between the plurality of clock signals to be measured . 

The clock skew estimator may include: a timing estimator 
15 operable to obtain a reference timing that is an edge timing of 
the reference signal and a tested timing that is an edge timing 
of each of the plurality of clock signals to be measured; a timing 
error estimator operable to obtain the timing difference between 
the tested timing and the reference timing; and a clock skew 
20 calculator operable to obtain the clock skew estimates between 
the plurality of clock signals to be measured from the timing 
difference obtained for each of the plurality of clock signals 
to be measured. 

25 The clock skew estimator may further include a corrector 

operable to correct the clock skew obtained by the clock skew 
calculator . 

The timing estimator may obtain a rising edge timing or a 
30 falling edge timing of eachof the reference signal and the plurality 
of clock signals to be measured. 

The timing estimator may include: an analytic signal 
transformer operable to transform each of the plurality of clock 
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signals to be measured into a complex analytic signal; an 
instantaneous phase estimator operable to obtain an instantaneous 
phase of the analytic signal; a linear instantaneous phase 
estimator operable to obtain a linear instantaneous phase of each 
5 of the plurality of clock signals to be measured based on the 
instantaneous phase obtained; and an initial phase estimator 
operable to obtain an ideal edge timing of each of the plurality 
of clock signals to be measured by obtaining an initial phase angle 
of the linear instantaneous phase. 



The timing estimator may further include: a linear trend 
remover operable to remove the linear instantaneous phase from 
the instantaneous phase to obtain an instantaneous phase noise; 
and a zero-crossing resampler operable to re-sample only data of 
15 the instantaneous phase noise around zero-crossing timings of a 
real part of the analytic signal and to output a timing jitter 
sequence of each of the plurality of clock signals to be measured. 



20 filter operable to receive each of the plurality of clock signals 
and to extract from the received clock signal frequency components 
around a fundamental frequency of the received clock signal thereby 
outputting a band-limited signal; and a Hilbert transformer 
operable to perform Hilbert transformation for the band- limited 

25 signal to generate Hilbert transform pairs of the clock signal 
to be measured. 

The analytic signal transformer may include: a time-domain 
to frequency-domain transformer operable to receive each of the 
30 clock signals to be measured and to transform the received clock 
signal into two-sided spectra in a frequency domain; a bandwidth 
limiter operable to extract from the two-sided spectra frequency 
components around a positive fundamental frequency thereof; and 
a frequency- domain to time-domain transformer operable to 
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The analytic signal transformer may include: a band-pass 
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inversely transform an output of the bandwidth limiter into a 
time-domain signal . 

The analytic signal transformer may include : a buffer memory, 
5 to which each of the plurality of clock signals to be measured 
is supplied, operable to store the supplied clock signal; an 
extracting portion operable to select and extract a section of 
the stored clock signal in such a manner that a section presently 
extracted partially overlaps a section previously extracted; a 

10 window function multiplier operable to multiply the extracted 
section by a window function; a transforming portion operable to 
transform the multiplied section into two-sided spectra in a 
frequency domain; a bandwidth limiter operable to extract, from 
the two-sided spectra transformed in the frequency domain, 

15 frequency components around a positive fundamental frequency of 
the supplied clock signal to be measured; an inverse transformer 
operable to inversely transforming an output of the bandwidth 
limiter into a time-domain signal; and an inverse window function 
multiplier operable to multiply the time-domain signal by an 

20 inverse of the window function to obtain the analytic signal that 
has been band-limited. 

The clock skew estimator may include ah analog- to-digital 
converter operable to receive the reference signal and each of 
25 the clock signals to be measured and to digitize the reference 
signal and each of the clock signals to be measured. 

The clock skew estimator may include a waveform clipper 
operable to receive the reference signal and each of the clock 
30 signals to be measured and to remove amplitude modulation 

components of the received clock signal to be measured to extract 
phase modulation components of the received clock signal. 

The analytic signal transformer may be operable to be 
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arranged to have a variable passband of each of the plurality of 
clock signals to be measured. 



The timing estimator may further include a low- frequency 
5 component remover operable to receive the instantaneous phase noise 
and to remove low frequency components of the instantaneous phase 
noise to output the instantaneous phase noise, from which the low 
frequency components have been removed, to the zero-crossing 
resampler . 

10 

According to the second aspect of the present invention, 
a clock skew measuring method for measuring a clock skew between 
a plurality of clock signals to be measured in a device under test, 
comprises : outputting the plurality of clock signals to be measured 

15 by selecting one of the plurality of clock signals to be measured 
one by one; and obtaining the clock skew estimates between the 
plurality of clock signals to be measured by measuring a timing 
difference between reference signal input to the device under test 
and each of the plurality of clock signals to be measured one by 

20 one . 

The reference signal may be a system clock signal supplied 
to the device under test . 

25 The outputting and selecting step may include determining 

which one of the plurality of clock signals is to be selected based 
on the reference signal. 

The clock skew estimating step may measure a deterministic 
30 component of the clock skew between the plurality of clock signals 
to be measured. 



The clock skew estimating step may measure a random component 
of the clock skew between the plurality of clock signals to be 



AD-0318 



.31 Q iuj ^ ■a.o-o m yia as: ^ ui .* 




measured. 



The clock skew estimating step may include: obtaining an 
edge timing of the reference signal as a reference timing; obtaining 
5 an edge timing of each of the plurality of clock signals to be 
measured as a tested timing; obtaining the timing difference 
between the tested timing and the reference timing; and obtaining 
the clock skew estimates between the plurality of clock signals 
to be measured from the timing difference obtained for each of 
10 the plurality of clock signals to be measured. 

The clock skew obtaining step may further include correcting 
the clock skew obtained from the timing difference. 

15 The obtaining of edge timing may obtain a rising edge timing 

or a falling edge timing of each of the reference signal and the 
plurality of clock signals to be measured. 

The timing estimating may include : transforming each of the 
20 plurality of clock signals to be measured into a complex analytic 
signal; obtaining an instantaneous phase of the analytic signal; 
obtaining a linear instantaneous phase of each of the plurality 
of clock signals to be measured based on the instantaneous phase 
obtained; and obtaining an ideal edge timingof eachof theplurality 
25 of clock signals to be measured by obtaining an initial phase angle 
of the linear instantaneous phase. 

The obtaining of edge timing may include : removing the linear 
instantaneous phase from the instantaneous phase to obtain an 
30 instantaneous phase noise; and re -sampling only data of the 

instantaneous phase noise around zero-crossing timings of a real 
part of the analytic signal to output a timing jitter sequence 
of each of the plurality of clock signals to be measured. 
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The transformation into the complex analytic signal may- 
include : extracting, from each of the plurality of clock signals 
to be measured, frequency components around a fundamental frequency 
of the clock signal to output a band-limited signal ; and performing 
Hilbert transformation for the band-limited signal to generate 
Hilbert transform pairs of the clock signal. 



The transformation into the complex analytic signal may 
include: transforming each of the plurality of clock signals to 
be measured into two-sided spectra in a frequency domain; 
extracting, from the two-sided spectra, frequency components 
around a positive fundamental frequency thereof; and inversely 
transforming the two-sided spectra that has been band- limited into 
a time-domain signal. 

The transformation into the complex analytic signal may 
include: storing eachof the clock signals tobe measured; selecting 
and extracting a section of the stored clock signal in such a manner 
that a section presently extracted partially overlaps a section 
previously extracted; multiplying the extracted sect ionby a window 
function; transforming the multiplied section into two-sided 
spectra in a frequency domain; extracting, from the two-sided 
spectra transformed in the frequency domain, frequency components 
around a positive fundamental frequency of the stored clock signal 
to be measured; inversely transforming the spectra that has been 
band-limited into a time-domain signal; and multiplying the 
time-domain signal by an inverse of the window function to obtain 
the analytic signal that has been band-limited. 

The obtaining of the timing difference between the tested 
timing and the reference timing may include: calculating a 
plurality of timing differences from the tested timing and the 
reference timing for each of the plurality of clock signals; and 
obtaining the mean value of the plurality of timing differences, 
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and wherein the clock skew obtaining obtains the clock skew between 
the plurality of clock signals to be measured based on the average 
of the plurality of timing differences. 



The clock skew obtaining may remove amplitude modulation 
components from the reference signal and each of the plurality 
of clock signals to be measured to extract phase modulation 
components thereof . 



The obtaining of the edge timing further may include removing 
low frequency components of the instantaneous phase noise. 

This summary of the invention does not necessarily describe 
all necessary features of the present invention so that the 
invention may also be a sub-combination of these described 
features. 



BRIEF DESCRIPTION OF THE DRAWINGS 

Fig . 1 illustrates an exemplary clock skew measurement using 
a time interval analyzer. 

Fig. 2A schematically shows a clock distribution network. 

Fig. 2B schematically shows timing analysis for clock skew. 

Fig. 3 shows an exemplary clock skew measurement circuit 
of the present invention. 

Fig . 4 shows another exemplary clock skewmeasurement circuit 
of the present invention. 

Fig. 5A schematically shows a clock distribution network. 

Fig. 5B schematically shows timing analysis for clock skew. 

Fig. 6A shows exemplary timing jitter Acp j [n] of a clock CLKj 
to be measured. 

5^P?Fig.\6B shows exemplary timing jitter Acp k [n] of a clock CLK k 
to be measured. 

Fig. 7 schematically shows a clock distribution network with 
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different clock domains. 

Fig. 8 schematically shows a principle of clock skew 
measurement using frequency multiplication by modulo M: M = 2. 

Fig. 9 shows an exemplary clock signal to be measured x(t) . 
5 Fig. 10 shows an exemplary analytic signal z(t). 

5j^2^> Fi^^ll shows an exemplary discontinuous instantaneous phase 

3^^/>F\g. 12 shows an unwrapped continuous instantaneous phase 

10 Fig. 13 shows an exemplary digitized clock signal x(t) . 

Fig. 14 shows an example of two-sided power spectra of the 
clock signal x(f) obtained by FFT. 

Fig. 15 shows an exemplary band-passed power spectra of the 
clock signal Z(f) . 
15 Fig. 16 shows an exemplary band-passed analytic signal z (t) 

obtained by inverse FFT. 

Fig. 17 shows an exemplary clock signal x(t) . 

Fig. 18 shows an exemplary analytic signal z (t) of the clock 
signal x (t) . 

20 Fig. 19 shows an exemplary instantaneous phase cp(t). 

Fig. 20 shows an exemplary instantaneous phase noise Acp (t) . 

U^21 sho^aja ^xem plary ^ t imin g- j i r i ^be y^ Afp [nl_ 
Fig. 22 shows exemplary adaptive zero-crossing points 
approximation . 

25 Fig. 23 shows an exemplary clock signal with AM components. 

Fig. 24 shows an exemplary clock signal without AM 
components . 

Fig. 2 5 shows an exemplary structure of a clock skew 
measurement apparatus of the present invention. 
30 Fig. 26 is a flowchart of an exemplary clock skew measurement 

method of the present invention. 

Fig. 2 7 is a flowchart of another exemplary clock skew 
measurement method of the present invention. 

Fig. 2 8 shows an exemplary structure of a timing estimator 



4> 



AD-0318 



11 



m OSS 2£ O S 



included in the clock skew measurement apparatus of the present 
invention. 

Fig. 2 9 is a flowchart of an exemplary timing estimation 
method included in the clock skew measurement method of the present 
5 invention. 

Fig. 3 0 shows an exemplary structure of an analytic signal 
transformer included in the clock skew measurement apparatus of 
the present invention. 

Fig. 31 is a flowchart of an exemplary analytic signal 
10 trans format ion method included in the clock skew measurement method 
of the present invention. 

Fig. 32 shows another exemplary structure of an analytic 
signal transformer included in the clock skewmeasurement apparatus 
of the present invention. 
15 Fig. 33 is a flowchart of another exemplary analytic signal 

transformation method included in the clock skewmeasurement method 
of the present invention. 

Fig. 34 shows still another exemplary structure of an 
analytic signal transformer included in the clock skewmeasurement 
20 apparatus of the present invention. 

Fig. 35 is a flowchart of still another exemplary analytic 
signal transformation method included in the clock skewmeasurement 
method of the present invention. 

Fig. 36 shows another exemplary structure of a clock skew 
25 measurement apparatus of the present invention. 

Fig. 3 7 is a flowchart of another exemplary clock skew 
measurement method of the present invention. 

Fig. 38 is a flowchart of still another exemplary clock skew 
measurement method of the present invention. 
30 Fig. 3 9 shows still another exemplary structure of a clock 

skew measurement apparatus of the present invention. 

Fig. 4 0 is a flowchart of further another exemplary clock 
skew measurement method of the present invention. 

Fig. 41 is a flowchart of further another exemplary clock 
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skew measurement method of the present invention. 

Fig. 42 shows another exemplary structure of a timing 

estimator included in the clock skew measurement apparatus of the 

present invention. 

Fig . 43 is a flowchart of another exemplary timing estimation 

method included in the clock skew measurement method of the present 

invention. 



system using the clock skew measurement apparatus of the present 
invention. 



The invention will now be described based on the preferred 
embodiments, which do not intend to limit the scope of the present 
invention, but exemplify the invention. All of the features and 
the combinations thereof described in the embodiment are not 
necessarily essential to the invention. 

Clock skew measuring method (1) 

First, clock skew is defined. With a clock source CLKg of 
a clock distribution network as a reference, as shown in Fig. 2A, 
the clock skew is given as a difference between T j cd and x k C d where 
i j cd and i k cd are delay times of distributed clock signals CLKj and 
CLKk required for reaching corresponding registers Rj and R k . 



Fig. 44 shows an exemplary arrangement of a clock skew testing 



DETAILED DESCRIPTION OF THE INVENTION 



1 Skew L cd L cd 



(1) 



Fig. 2B shows a timing of the clock skew. 



The delay times x j C d and i k C d are respectively represented 



as follows: 



L cd l cd L cd 



(2) 
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(3) 



where rising edge timings of the clock signals CLKg, CLKj and CLK k 
are t^cci, t^cd. 

and t k cd/ respectively. As is apparent from those, 
5 the clock skew between the distributed clock signals CLKj and 

CLKk can be obtained from the time difference between the edge 
timing of each of the distributed clock signals CLKj and CLK k and 
the reference timing, where the edge timing of the clock source 
CLKg is the reference timing. 



10 
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the clock skew measuring method of the present embodiment , 
the distributed clocks to be measured are selected and brought 
out to the outside of the chip by means of a clock signal selector, 
for examfole, a multiplexer. Then, the time difference between 
15 the edge timing of the clock signal and the reference timing is 
measured fbr each of the selected clock signals, so that an error 
of the time differences thus measured is obtained, thereby 
measuring the clock skew. For simplification, the method for 
measuring thet skew between two distributed clock signals is 
20 described belc 

Fig. 3 shows an exemplary clock skew measuring circuit of 
the present embodiment. The clock skew measuring circuit includes 
buffers for selecting and outputting the distributed clocks CLKj 

25 and CLK k to an output pin, and a multiplexer for selecting the 
distributed clock. The buffer can select whether or not the 
associated distributed clock is to be output to its output in 
accordance with an external input signal (ENBsignal) . Thebuffers 
BUFj and BUF k have identical design, and the transmission delay 

30 time in each buffer is assumed to be d BU F- In addition, it is assumed 
that delay times of lines Pj and P k , the transmission delay time 
of the multiplexer, and the delay time from the clock signal 
selecting element to the output are dj , d k , dnux and d OU r/ 
respectively . 
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The clock skew measuring method of the present embodiment 
uses the ENB signal that is synchronized with the edges of the 
clock source CLKg of the distributed clocks as the reference signal . 
5 Since the ENB signal is in synchronization with CLKg, a time 
difference between an edge timing t e c n d b of the ENB signal and the 
rising edge timing t s cd of CLKg is always constant. 

l cd L cd u const v ' 

10 

Next, a procedure of the clock skew measuring method of the 
present embodiment is described. First, the clock skew measuring 
method of the present embodiment sets a selection signal SEL for 
the multiplexer to zero to select CLKj , thereby measuring a 
15 difference Dj between the edge timing of the ENB signal and that 
of CLKj output to the output pin. Dj is obtained as follows: 

Dj={tL+d BUF +dj+d MUX +d OUT yC (5) 

20 where the rising edge timings of the ENB signal and CLKj are t^f 
and t J cd , respectively. 

Next, the selection signal SEL for the multiplexer is set 
to one so that CLK k is selected, and a difference D k between the 
25 edge timing of the ENB signal and that of CLK k output to the output 
pin. D k is obtained as follows: 

D k = (t k cd + d BUF + d h + d MVX + d OUT )- C ( 6 ) 

30 where the rising edge timings of the ENB signal and CLK k are t% b 
and t k cd , respectively. 



1003318B „05S30e 

14 




AD-0318 



Finally, a difference between the thus measured differences 
D k and Dj . 

D k - D. = fc + d BUF + d k + d MUX + d OUT )- C }~ {fa + d BUF + dj + d MUX + d 0UT )-C } 

5 

When this equation is calculated by using Equation (4), 

10 is obtained. Using Equations (1) , (2) and (3) , the difference 
between D k and Dj is obtained as follows. 

D k -Dj=(r k cd -ri d )-(d k -dj)=T£, -{d k -dj) (7) 

15 Therefore, when the lines Pj and P k that respectively connect 

the buffers to the multiplexer are laid out in such a manner that 
the delay times dj and d k are equal to each other, the clock skew 
between CLKj and CLK k can be obtained by obtaining the difference 
between D k and Dj by using the clock skew measuring method of the 

20 present invention. 

Moreover, if the delay times dj and d k are not equal to each 
other, the clock skew between CLKj and CLK k can be obtained by 
correcting the difference between D k and Dj . 

25 

Tj* w =(D k -D J )+(d k -d J ) (8) 

In the above, the difference between d k and dj may be obtained 
by circuit simulation or actual measurement using, for example, 
30 a circuit shown in Fig. 4. 

Fig. 4 shows a modification of the clock skew measuring 
circuit of the present embodiment. In Fig. 4, propagation delay 
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times D enb and Df nb that are times required for propagating the ENB 

J K 

signal until the output pin OUT via the lines Pj andP k , respectively, 
can be measured by supplying the ENB signal while the MEAS signal 
is set to one. 

5 

D ent = + deuF + d _ + d ^ + dofjT y t enb = ^ + ^ + ^ + ^ 

at 6 = (c + ^ buf +d k + d MUX + d OUT )-c = ^ fi£/F + d k + d MUX + d OUT 

The difference between Z>" 6 and D e k nb is obtained as follows. 

10 

DT b -Df=d k -dj 

Thus, the difference between dj and d k can be obtained by obtaining 
the transmission times D e " b and D e k nb that are the times required 
15 for the transmission of the ENB signal to the output pin OUT via 
the lines Pj and P k by using the circuit shown in Fig. 4 # and then 
calculating the difference between the thus obtained propagation 

delay times Z^and D e k nb . 

20 Clock skew measuring method (2) 

Next , a clock skewmeasuringmethodusing instantaneous phase 
estimation by an analytic signal is described. This method can 
obtain a deterministic component and a random component of the 
clock skew. 

25 

First , the clock skew between the clock signals having j itter 
is defined. When the clock source CLKg of the clock distribution 
network is assumed to be the reference, for example, as shown in 
Fig. 5A, the clock skew is given as the difference between the 
30 delay times i j C d and i k cd that are the times required for the 

propagation delay of the distributed clock signals CLKj and CLK k 
to reach the corresponding registers Rj and R k . Since the rising 
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edge timing of each clock signal is varied for each cycle due to 
the jitter, the clock skew of each cycle, T^{nT) , is expressed 
as follows. 

TiL{nT)=T k cd {nT)-ri d {nT) (9) 

Fig. 5B shows the timings of the clock skew. In this description 
and Fig. 5B, T is a fundamental period of the clock signal to be 
measured. 



The rising edge timings of the clock signals CLKg, CLKj and 
CLK k are assumed to be t 9 C d, t j C d and t k cd/ respectively. When ideal 
clock edge timings of the clock signals CLKg, CLKj and CLK k (i.e. , 
the clock edge timings of the clock signals with no jitter) are 
15. assumed to be (nT) g , (nT)j and (nT) k/ respectively, the delay times 
of each cycle i j cd (nT) and i k cd (nT) are respectively expressed as 
follows . 

^(«r)=^(«r)-^(«r) 

20 = [tU (nT)- (nT)j]-[t? d (nT)- (nT) J+ {(nT)j - (nT) g } 



- t s.j . 

'■delay T 



a ^^^)~ a ^M^)] tsec] (io) 



T k cAnT) = t k cd {nT)-t* cd {nT) 

= iL inT)-{nT\\- [,* (nT)-(n T \ ] + {{nT) k - (n T \ } 



25 =r a % + 



[sec] (11) 



In the above expressions, 

r g J ay ={nT) j -{nT) g [sec] (12) 
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(13) 



are respectively time differences between the ideal clock edge 
timings of the distributed clocks CLKj and CLK k and the ideal clock 
5 edge timing of the clock source CLKg, and correspond to the 
deterministic component s of the propagation delays that are defined 
by transmission paths (deterministic propagation delay times) . 
Moreover, Acp 9 [n] (T g /2n. (=t g cd (nt) - (nT) g ) , Acp j [n] (Tj/2n. (=t j cd (nt) 
- (nT)j) andAcp k [n] (T k /2n. (=t k cd (nt) - (nT) k ) respectively represent 
10 timing jitter sequence of the clocks CLKg, CLKj and CLK k (unit is 
second) . How to estimate the timing jitter of the clock signal 
will be described later. From Equations (9) , (10) and (11) , the 
clock skew Tgjn] ( = T£^(nT)) between CLKj and CLK k is estimated as 
follows . ■ 

15 




[sec] (14) 



20 The first term of Equation (14) 

* J sL=(»T) k -(nT)j [sec] (15) 

25 represents the difference of the rising edge timing between the 
ideal clocks of CLKj and CLK k/ and is an example of the deterministic 
component of the clock skew determined from the paths in the clock 
distribution network, for example. In addition, the second term 
of Equation (14) 
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is an example of the random component of the clock skew due to 
5 the timing jitter contained in the respective clock signals. 

The deterministic clock skew value can be estimated by- 

obtaining, for example, the instantaneous phases of two signals 
CLKj and CLK k and then obtaining the difference between linear 
10 phase components of the instantaneous phases of these two signals . 
The fundamental -frequency cosine components of CLKj and CLK k are 
assumed as follows. 

Xj {f) - Aj cos(^ J (t)) = Aj cos 
15 x k (t) = A k cos(^* A k cos 

Here, the instantaneous phases of Xj (t) and x k (t) are represented 
as the sum of the linear instantaneous phase component 2nt/T L 
containing the fundamental period T L (L = j , k) , the initial phase 

20 angle (L = j , k) , and the instantaneous phase noise component 
A0 L (t) (L = j , k) . 

</> k (t) = ^-t + ti-A0 k {t) [rad] (19) 

25 

The estimation method of the instantaneous phase of the clock signal 
will be described later . When A^(/) = 0 in Equations (18) and (19), 



2n_ 
T 



t + M-Wit) 



\ 2 k J 



(16) 



(17) 
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the linear instantaneous phases of the clocks with no jitter 



*LrM = -fr' + ^ f rad 3 (20) 



are obtained. In this case, the ideal rising edge timings of CLKj 
andCLK k , t= (nT) j , (nT) k/ are times at which the linear instantaneous 
phase components in left -hand side of Equations (18) and (19) become 
(2nn - n/2) , respectively, and therefore the following 
10 relationships are satisfied. 

{nT). = [ 2w ^~f ~^o)-^ [sec] (22) 



\ T 

[sec] (23) 

)2n 



Thus, from Equation (15) , the deterministic clock skew value 




is obtained. In general, the fundamental periods of the 
distributed clock signals CLKj and CLK k are equal to each other 
(Tj = T k ) . That is, the deterministic clock skew value between 
25 the two signals to be measured can be obtained as the difference 
of the initial phase angle in the linear instantaneous phase between 
the two signals to be measured. 
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The initial phase angle 0 O can be obtained by performing 
a least -square fit of a to the instantaneous phase data <p(k) and 
then determining 0 O that makes 



(25) 



be minimum. At this time, the initial phase angle is given as 
follows . 



10 



2(2^ + l)2#(*)-6X*^(*) 



N(N-\) 



k=\ 



(26) 



15 



Moreover, the initial phase angle </> 0 of the signal to be 
measured x(t) can be obtained by performing cosine wave fitting 
using the least -squares method for the clock waveform data x(k) 
or its fundamental sine component and then estimating 0 O that makes 



Y\x{k)-Aco$\ ^-k + 0 o 



k=\ ' 



V T 



(27) 



be maximum in accordance with a maximum likelihood estimation 
20 method. In this case, the initial phase angle is given by the 
following expression . 



= -arctan 



f N 2tt ^ 

]^;c(£)sin k 

k=\ T 

^jc(£)cos — k 

\k=\ T 



(28) 



25 In the above, it is assumed that the corresponding clock 

edges of the two signals to be measured are not separated by one 
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period or more . On the other hand , in a case where the corresponding 
clock edges are away from each other by one period or more, the 
deterministic clock skew value can be obtained as the sum of the 
difference of the initial phase angle and the offset time of the 
5 clock edges. 



10 source has a strong relationship with the clock signal of the clock 
source. As a result of this, the phase noise of the distributed 
clock signal (timing jitter sequence) generally has a similar 
tendency to that of the phase noise of the clock source (timing 
jitter sequence). Thus, the timing jitter sequence of the 

15 distributed clock signals obtained by distribution from the same 
clock source have the same tendency (see Figs. 6A and 6B) . 
Therefore , the offset amount n of f set of the corresponding clock edges 
of the two signals to be measured can be estimated by obtaining 
correlation between the timing jitter sequence thereof and then 

20 finding a value of the offset that provides the maximum correlation . 
Also, the offset amount n 0 ff ge t may be obtained from the value of 
offset that provides the maximum correlation of the instantaneous 
phase noises. 

25 Moreover, the deterministic clock skew value maybe obtained 

by obtaining zero-crossing times of the signals to be measured 
and then calculating an average value of the time differences 
between the corresponding zero-crossing times. 

30 Next, the clock skew between the clock signals having 

different frequencies is described. In this description, the 
clock distribution network shown in Fig. 7 is considered. A clock 
source PLL g multiplies a system clock CLKg that is supplied from 




[sec] 



(29) 



The clock signal obtained by distribution from the clock 
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the outside, by M, and then distributes the clocks CLKj and CLK k 
to the network. In Fig. 8, (a) shows the system clock CLKg while 

(c) shows the clock CLKj obtained by multiplication. A9[l] [rad] 
of the system clock CLKg represents timing fluctuation with respect 
5 to the ideal clock edge of that edge. Therefore, when edges of 
the ideal clock obtained by multiplication by M are considered, 
as shown in Fig. 8(b) , and then A9[l] is copied to obtain (M-l) 
copies, A©[[/z/mJ] corresponds to A^ y [«] one by one. Please note 
that [xj represents largest integer less than or equal to x. When 
10 the clock skew between CLKj and CLKg is obtained by using Equation 

(14) , Equation (30) is obtained. 




-A© 



n 

~M 



±G_ 

In 



[sec] 



(30) 



15 



20 



The deterministic clock skew value between CLKj and CLKg 

is represented as the time difference between the ideal clock edge 
of CLKj, (nMT) j, and the ideal clock edge of the system clock CLKg, 
( nMT ) g / and can be obtained from the initial phase angles of the 
respective clocks by the following equation. 



T C s iL=(nMT) j -{nMT) c 



f T.\ (MT 0 



[sec] 



(31) 



Since the clock CLKj is obtained by multiplication of the system 
25 clock CLKg by M, the fundamental period of CLKg, T g , is equal to 
M times the fundamental period Tj of CLKj (T G = MTj) . 

Next, a procedure for obtaining the clock skew between the 
distributed clocks CLKj and CLKg by using the reference clock signal 
30 CLK R is described. 
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First, only CLKj and CLK R are sampled simultaneously, and 
the skew between CLKj and CLK R/ 



„H K^ii 
JL ixJl ^ -^U 



24 






[sec] 



(32) 



is obtained by using Equation (14) . Then, only CLK k and CLK R are 
simultaneously sampled, and the skew between CLK k and CLK R , 



10 



15 



T R,k r ] _ , 

1 Skew l n J ~~ c Skew ^ 



[sec] 



(33) 



is obtained. Finally, by obtaining the difference between the 
clock skew sequence obtained as mentioned above, the clock skew 
between CLKj and CLK k is obtained as follows. 



1 L Skew T 



i 6 Stew ^ 



=(^-^) + {h'["6)- A ^K3-h /[ "&)-^' [ "(^)] 



[sec] 



(34) 



20 The above procedure can be applied to the clock signals having 

different frequencies. Thus, the clock skew measuring method of 
the present embodiment can measure the clock skew between the 
on-chip clock signals by using the system clock supplied to the 
semiconductor integrated circuit under test from the outside 

25 thereof as the reference clock signal. 

The clock skew measuring method of the present embodiment 
performs simultaneous sampling for the distributed clock CLKj and 
the system clock CLKg and then obtains the skew between CLKj and 
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5 by using Equation (30) . Next, the other distributed clock CLK k 
and the system clock CLKq are simultaneously sampled, and the skew 
between CLK k and CLKg, 



10 



rG,k 
l Skew 
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is then obtained similarly. Finally, by calculating the 
difference between the clock skew sequence obtained as mentioned 
above, the clock skew between CLKj and CLK k , 



15 T^[n] = T™[n]-T^[n] 



I G,k 
I 1 Skew T 
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20 



[sec] (37) 



is estimated. In addition, as described in the description of 
the first clock skew measuring method, when the delay times dj 
and d k of the clock-output lines are not equal to each other, the 
25 clock skew between CLKj and CLK k can be obtained with high precision 
by obtaining the difference between dj and d k by simulation or 
the like and then correcting the clock skew obtained by Equation 
(37) by using the obtained difference between dj and d k . 
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[sec] (38) 

5 

As a result, according to the clock skew measuring method 
of the present embodiment, the measurement can be performed by 
selecting and making the on-chip clock signals distributed over 
the semiconductor chip output to the outside of the chip one by 

10 one. Thus, by adding the clock signal selecting element for 
selecting the clock signal and outputting it to the outside to 
the device under test, the number of high-frequency clock- output 
pins, that may increase the cost, can be greatly reduced. In the 
present embodiment, for example, the number of high-frequency 

15 clock-output pins can be reduced to be minimum. Therefore, the 
clock skew measuring method is suitable for analysis or test for 
VLSIs. 

The clock skew measuring method of the present embodiment 
20 can estimate not only the clock skews between the distributed clock 
signals of MPU as described above but also the clock skews between 
other types of signals. 

In the above description, a case where the output of the 
25 clock signal selecting element is one is described as an example. 
However, the present invention is not limited to the above case. 
The present invention can be applied to the clock skew measurement 
using the clock signal selecting element having two or more outputs . 

30 Instantaneous phase estimation using analytic signal 

An analytic signal z(t) of a real signal x(t) is defined 
by the following complex signal. 
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(39) 



In the above, j is an imaginary unit and the imaginary part x(t) 
5 of the complex signal z(t) is Hilbert transform of the real part 
x(t) . 



10 



Hilbert transform of a time-variant waveformx (t) is defined 
as follows. 



*(0=#[*(0]=l7r^ 

71 J J-T 



(40) 



In the above, x(t) is convolution of functions x(t) and (1/nf) . 
That is, Hilbert transform is equivalent to the output obtained 
15 whenx(t) is made to pass through a full band-pass filter . However, 
the output x(t) has the spectra component but the phase shifted 
by n/2 . 

The instantaneous phase waveform (f>(f) of the real signal x ( t ) 
20 can be obtained from the analytic signal z(t) by using the following 
equation. 



<p(t)= tan" 1 



(41) 



25 Next , algorithm for estimating the instantaneous phase using 

Hilbert transformation is described. First, the signal to be 
measured shown in Fig. 9, 



30 



At) -A 



f 

cos 



2?r 



(42) 
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is subjected to Hilbert transformation so as to obtain a signal 
corresponding to the imaginary part of the complex signal, 
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= = ^sin —t + <f> 0 - &<p(t) 



(43) 



Thus, the signal to be measured, x(t) , is transformed to the analytic 
signal z(t) as follows. 



10 



15 



z(f) = jc(V)+ = y4 cos — / + (f> 0 - k<p(t) I + jA sin 



2^ 



~t + 



<f> Q -&(f>{t) 



(44) 



The analytic signal obtained by transformation is shown in Fig. 
10 . The obtained analytic signal has been subjected to band-pass 
filtering. This is because the jitter corresponds to the 
fluctuation of the fundamental frequency of the signal to be 
measured and therefore only the signal components around the 
fundamental frequency of the signal to be measured are used in 
jitter analysis. Then, a phase function (/){t) is estimated from 
the obtained analytic signal z (t) by using Equation (41). 



20 



#(0= 



mod27r 



[rad] 



(45) 



The phase function <f>{i) is represented using a principal value of 
the phase in the range from -n to +n and has a discontinuous point 
in the vicinity of a point at which the phase changes from -n to 
25 +n. The estimated phase function </>{f) isshowninFig.il. Finally, 
by unwrapping the discontinuous phase function <j>(t) (that is, 
appropriately adding an integral multiplication of 2n to the 
principal value <j){t)) , the discontinuity can be removed, thereby 
the continuous instantaneous phase (f>{t) can be obtained. 

30 
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f(t) = ±?Lt + <t> 0 - A#(t) [ rad] (46) 

The continuous instantaneous phase function </)(t) after being 
unwrapped is shown in Fig. 12. 

5 

Transformation into analytic signal using Fast Fourier 
Transformation 

The transformation from the real signal to the analytic 
signal can be achieved by digital signal processing using Fast 
10 Fourier Transformation (FFT) . 

First , FFT is applied to the digitized signal to be measured, 
x(t) , shown inFig. 13, so that two-sided spectra (which has positive 
and negative frequencies) x(f) of the signal to be measured is 

15 obtained. The obtained two-sided spectra x(f) is shown in Fig. 
14. Then, only data around the fundamental frequency in the 
positive frequency components of x(f ) is made to remain while all 
the other data is replaced with zero, and the positive frequency 
components are doubled. These operations in the frequency domain 

20 correspond to band limiting for the signal and transformation to 
the analytic signal in the time domain. The thus obtained signal 
in the frequency domain, Z(f) , is shown in Fig. 15. Finally, the 
obtained signal Z(f) is subjected to inverse FFT, thereby the 
analytic signal z(t) that has been band limited can be obtained. 

25 The analytic signal after being band limited, z(t) , is shown in 
Fig. 16. 

Moreover, in a case where the object is to estimate the 
instantaneous phase, the operation for doubling the positive 
30 frequency components can be omitted. 

Timing jitter estimation 

Next, a method for estimating the timing jitter, that is 
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used in the clock skew measuring method of the present embodiment, 
is described. 



The clock signal having no jitter is a square wave having 
5 a fundamental frequency f 0 . This signal can be decomposed into 
harmonics composed of frequencies of f 0 , 3f 0 , 5f 0 , . - . by Fourier 
analysis. Since the jitter corresponds to the fluctuation of the 
fundamental frequency of the signal to be measured, only the signal 
components around the fundamental frequency are considered in the 
10 jitter analysis. 

The fundamental -frequency sinusoidal wave component of the 
clock signal having jitter (signal to be measured) is represented 
as follows: 

15 

x(t) = A cos(0(* )) = A cos 

where the amplitude and the fundamental period are A and T 0 , 
respectively. In the above, (f>(t) is the instantaneous phase of 
20 the signal to be measured and is represented by the sum of the 
linear instantaneous phase component 2nt/T 0 containing the 
fundamental period T 0 , the initial phase angle cp (which can be 
made to be zero in calculation) and the instantaneous phase noise 
component Acp(t) . 

25 

When the instantaneous phase noise component Acp (t) is zero, 
the zero-crossing points at which the signal to be measured rises 
are away from each other with constant periods T 0 - Acp(t) that 
is not zero causes the zero-crossing point of the signal to be 
30 measured to fluctuate. That is, Acp(nT 0 ) at the zero-crossing point 
nT 0 represents the time-variant fluctuation at the zero-crossing 
point, and is called as timing jitter. Thus, the timing jitter 
of the signal to be measured can be obtained by estimating the 
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instantaneous phase cp(t) of the signal to be measured and then 
obtaining the difference between the instantaneous phase at the 
zero-crossing point and the linear phase (that corresponds to the 
phase waveform of the ideal clock signal having no jitter) , 2nt/T 0 
5 + cp, that is, the instantaneous phase noise Acp(t). 



According to the timing jitter estimation method of the 
present embodiment, the signal to be measured shown in Fig. 17 
is first transformed into the analytic signal z (t) . The analytic 

10 signal z(t) after transformation is shown in Fig. 18. In Fig. 
18, the solid line represents the real part of the analytic signal 
while the broken line represents the imaginary part thereof . Next , 
the instantaneous phase cp(t) of the signal to be measured is 
estimated from the analytic signal z (t) . The instantaneous phase 

15 waveform cp (t) obtained by the estimation is shown in Fig . 19. Next, 
a least squares fit of a straight line is performed for the 
instantaneous phase data, thereby obtaining the linear 
instantaneous phase cpunear (t ) that corresponds to the instantaneous 
phase waveform of the ideal signal having no jitter. Then, the 

20 instantaneous phase noise component of the signal to be measured, 
Acp(t), is obtained by calculating the difference between the 
instantaneous phase cp(t) and the linear instantaneous phase 
<Piinear(t) . The obtained instantaneous phase noise waveform Acp(t) 
is shown in Fig. 20 . Then, the instantaneous phase noise waveform 

25 Acp(t) is sampled at timings (approximated zero-crossing points) 
which are the closest timings to the respective zero-crossing 
points of the real part x(t) of the analytic signal z(t) , and the 
instantaneous phase noise at the zero-crossing timing nT 0 , that 
is, the timing jitter Acp [n] (= Acp(nT 0 )) is then estimated. The 

30 waveform of the estimated timing jitter, Acp [n] , is shown in Fig. 
21 . 

In the timing jitter estimation method of the present 
embodiment, the timing jitter may be estimated by removing 
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amplitude modulation (AM) components of the signal to be measured 
by means of a waveform clipper so as to make only phase modulation 
(PM) components corresponding to the jitter remain, thereby the 
high-precision jitter estimation can be performed. 



Moreover, in the timing jitter estimation method of the 
present embodiment, low- frequency components of the phase noise 
signal maybe removed by means of a low-frequency component remover. 

10 Detection method of approximated zero-crossing point 

Next, how to detect the approximated zero-crossing point 
is described. First, assuming that the maximum value and the 
minimum value of the real part x(t) of the analytic signal of the 
input signal to be measured are values at 100% level and 0% level, 

15 respectively, a signal value at 50% level, V 50 %, is calculated as 
a level of the zero cross. Then, differences between adjacent 
sampled values of x (t) and the 50%-level value V 50 %/ (x(j-l) - V 50 %) 
and (x(j) - V 50 %) , are obtained and thereafter the product of these 
differences, (x(j-l) - V 50 %) x (x(j) - V 50 %) is calculated. When 

20 x(t) goes across 50% level, i.e., the zero-crossing level, signs 
of these sampled values (x(j-l) - V 50 %) and (x(j) - V 50 %) change 
from minus to plus or from plus to minus. Thus, when the above 
product is negative, it is determined that x(t) goes across the 
zero-crossing level . Therefore, one of the times (j-1) and j which 

25 corresponds with one of the sampled values (x(j-l) - V 50 %) and (x(j) 
- V 50 %) that has a smaller absolute value than the other at that 
time, is obtained as the approximated zero-crossing point. Fig. 
22 shows the waveform of the real part x(t) of the analytic signal. 
A mark in Fig. 22 represents the closest point (approximated 

30 zero-crossing point) to the rising zero-crossing point detected. 

Waveform clipping 



input thereto and makes only the PM components corresponding to 



5 



A waveform clipper removes the AM components from the signal 



AD-0318 



.JL U « JSi. JUHi O uu SO =» = 1UH e 

the jitter remain. The waveform clipping is achieved by 1) 
multiplying the input signal which is analog or digital by a 
constant; 2) replacing a signal value that is larger than a 
predetermined threshold value 1 with the threshold value 1; and 
5 3) replacing a signal value that is smaller than another 

predetermined threshold value 2 with the other threshold value 
2. Here it is assumed that the threshold value 1 is larger than 
the threshold value 2 . The clock signal containing the AM 
components is shown in Fig. 23. Since an envelope of the 
10 time-variant waveform is varied, the AM components are understood 
to exist. On the other hand, Fig. 24 shows the clock signal after 
being clipped by the waveform clipper. The time-variant waveform 
in Fig. 24 has a constant envelope and therefore it is confirmed 
that the AM components are removed. 

15 

Hereinafter, another example of the clock skew measuring 
apparatus of the present embodiment is described. For 
simplification, an example of the apparatus and method for 
measuring the clock skew between two signals to be tested is 
20 described below. However, the present invention can measure the 
clock skews between three or more clock signals in a similar manner . 




Fig . 25 illustrates another exemplary structure of the clock 
skew measuring apparatus of the present embodiment . The clock 

25 skew measuring apparatus 100 includes buffers 101 having signal 
inputs ENB for selecting whether or not a plurality of clock signals 
to be tested in a device under test are to be output to the outside 
of the chip, respectively, a clock signal selector 102 for 
outputting the clock signals to be tested selecting one of the 

30 clock signals input thereto one by one, which has a clock signal 
selecting input SEL; and a clock skew estimator 103 for measuring 
a timing difference between a reference signal input to the device 
under test and the clock signal selected by the clock signal selector 
for each of the clock signals to be tested, thereby obtaining the 
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clock skew between the clock signals to be tested. The clock skew 
estimator 103 includes a timing estimator 104a which obtains edge 
timing of the reference signal (reference edge timing) , another 
timing estimator 104b which obtains edge timing of the clock signal 
5 to be tested (tested timing) , a timing error estimator 105 which 
obtains the timing difference between the tested timing and the 
reference timing, a clock skew calculator 106 which obtains the 
clock skew between the clock signals to be tested from the timing 
difference obtained for each of the clock signals to be tested, 
10 and a corrector 107 which corrects a value of the clock skew obtained 
by the clock skew calculator 106. A specific structure of the 
timing estimator will be described later. 



15 measuring apparatus 100 of the present embodiment is described. 
Fig. 26 shows an exemplary procedure of the clock skew measuring 
method of the present embodiment. First, in Step 2 01, "0" is 
supplied to the selecting input SEL of the clock signal selector 
102, thereby selecting the clock signal CLK0 . Then, in Step 202, 

20 a signal input which enables the buffer output. is supplied to the 
ENB inputs of the buffers 101 as a reference signal of the clock 
skew measurement in synchronization with the reference clock edge 
of the device under test . Then, the timing estimator 104a measures 
the edge timing of the reference signal ENB in Step 203, and the 

25 timing estimator 104b measures the edge timing of the clock signal 
CLKO, that is output to the clock output pin of the device under 
test, in Step 204. In Step 205, the timing error estimator 105 
obtains the timing difference between the tested timing measured 
in Step 204 and the reference timing measured in Step 203. Next, 

30 in Step 206, "1" is supplied to the selecting input SEL of the 
clock signal selector 102 , thereby selecting the clock signal CLK1 . 
Then, in Step 2 07 , the signal input which enables the buffer output 
is supplied to the ENB input of the buffers 101 as the reference 
signal of the clock skew measurement in synchronization with the 



Next, a clock skew measuring operation of the clock skew 



AD-0318 




reference clock edge of the device under test. Then, the timing 
estimator 104a measures the edge timing of the reference signal 
ENB in Step 2 08, and the timing estimator 104b measures the edge 
timing of the clock signal CLK1, that is output to the clock output 
5 pin of the device under test , in Step 209. In Step 210, the timing 
error estimator 105 obtains the timing difference between the 
tested timing measured in Step 2 0 9 and the reference t iming measured 
in Step 208. Next, in Step 211, the clock skew calculator 106 
calculates the clock skew between CLK0 and CLK1 by obtaining the 

10 difference between the timing differences measured in Steps 2 05 
and 210. Finally, in Step 212, the corrector 107 corrects the 
clock skew obtained in Step 212, thereby finishing the procedure. 
In Steps 205 and 210 in which the difference between the tested 
timing and the reference timing is obtained, the timing error 

15 estimator 105 obtains that timing difference in accordance with 
Equations (5) and (6) . Moreover, in Step 211 where the clock skew 
between CLK0 and CLK1 is obtained, the clock skew calculator 106 
calculates the difference between the timing differences by using. 
Equation (7) . In addition, in Step 212 for correcting the obtained 

20 clock skew, the corrector 107 corrects the clock skew with the 
difference between the delay times of the lines (wires) by using 
Equation (8) . In step 211, the clock skew calculator 106 may obtain 
the absolute value of Equation (7) , if necessary. Moreover, Step 
212 may be omitted when the lines for outputting the clocks are 

25 designed and layout in such a manner that the difference between 
the delay times of the lines is zero. Furthermore, in order to 
improve the precision of the clock skew measurement, Steps 201 
to 212 may be repeated a plurality of times so that an average 
of the obtained clock skews can be obtained. 

30 

Fig. 2 7 shows another exemplary procedure of the clock skew 
measuring method of the present invention. First, in Step 301, 
"0" is supplied to the selecting input SEL of the clock signal 
selector 102, thereby selecting CLK0 . Then, inStep302, the timing 
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estimator 104a measures the system clock input to the device under 
test, so as to obtain the reference edge timing. In Step 303, 
the timing estimator 104b measures the edge timing of the clock 
signal CLKO that is output to the clock output pin of the device 
5 under test . The timing error estimator 105 then obtains the timing 
difference between the tested timing measured in Step 303 and the 
reference timing measured in Step 3 02 , in Step 304. Next, in Step 
305, "1" is supplied to the selecting input SEL of the clock signal 
selector 102, thereby selecting CLK1 . Then, inStep306, the timing 

10 estimator 104a measures the system clock so as to obtain the 
reference edge timing. In Step 307, the timing estimator 104b 
measures the edge timing of the clock signal CLK1 that is output 
to the clock output pin of the device under test . The timing error 
estimator 105 then obtains the timing difference between the tested 

15 timing measured in Step 3 07 and the reference timing measured in 
Step 306, in Step 308. Then, in Step 309, the clock skew calculator 
106 obtains the clock skew between CLKO and CLK1 by obtaining the 
difference between the timing differences measured in Steps 304 
and 308. Finally, in Step 310, the corrector 107 corrects the 

20 clock skew obtained in Step 3 09, thereby finishing the procedure. 
In Steps 304 and 308 where the timing difference between the tested 
timing and the reference timing is obtained, the timing error 
estimator 105 obtains that timing difference, that is, the skew 
between the system clock and the clock signal to be measured, in 

25 accordance with Equations (3 5) and (36) . Moreover, in Step 309 
where the clock skew between CLKO and CLK1 is obtained, the clock 
skew calculator 106 calculates the difference between the timing 
differences by using Equation (37) . In addition, in Step 310 for 
correcting the obtained clock skew, the corrector 107 corrects 

30 the clock skew with the difference between the delay times of the 
lines by using Equation (38) . In step 309, the clock skew 
calculator 106 may obtain the absolute value of Equation (37) , 
if necessary. Moreover, Step 310 may be omitted when the lines 
for outputting the clocks are designed and laid out in such a manner 
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that the difference between the delay times of the lines is zero. 
Furthermore, Steps 302, 303, 306 and 307 for obtaining the edge 
timings of the clock signals may be replaced with a procedure shown 
in Fig. 29. 



The timing estimators 104a and 104b mentioned above can be 
arranged as shown in Fig. 2 8 . Fig. 2 8 shows an exemplary structure 
of the timing estimator of the present embodiment. The timing 
estimator 400 includes an analytic signal transformer 401 which 

10 transforms the clock signal input thereto into a complex analytic 
signal; an instantaneous phase estimator 402 for obtaining an 
instantaneous phase of the analytic signal ; a linear instantaneous 
phase estimator 403 for obtaining a linear instantaneous phase 
of the input signal from the instantaneous phase; an initial phase 

15 estimator 404 for obtaining an ideal edge timing of the input clock 
signal by obtaining an initial phase angle of the linear 
instantaneous phase; a linear trend remover 4 05 for removing the 
linear instantaneous phase from the instantaneous phase so as to 
obtain an instantaneous phase noise ; and a zero-crossing resampler 

20 406 for re-sampling only the instantaneous phase noise data around 
zero-crossing timings of the real part of the analytic signal and 
outputting the timing jitter sequence of the input signal. The 
analytic signal transformer 401 can have the structure shown in 
any of Figs. 30, 32 and 34. 



Next , an operation of the timing estimator 400 of the present 
embodiment for estimating the edge timing of the clock signal to 
be measured is described. Fig. 2 9 shows an exemplary procedure 
of the timing estimation in the present embodiment. First, in 
30 Step 501, the analytic signal transformer 401 transforms the input 
clock signal to be measured into the analytic signal containing 
predetermined frequency components allowed to pass therethrough . 
Then, in Step 502 , the instantaneous phase estimator 402 estimates 
the instantaneous phase of the signal to be measured by using the 
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analytic signal obtained by the analytic signal transformer 401. 
The linear phase estimator 403 then estimates the linear 
instantaneous phase corresponding to an ideal clock signal from 
the instantaneous phase estimated by the instantaneous phase 
5 estimator 402, in Step 503. Then, in Step 504, the initial phase 
estimator 404 obtains the ideal edge timing of the input clock 
signal by obtaining the initial phase angle of the linear 
instantaneous phase estimated by the linear phase estimator 403. 
In Step 505, the linear trend remover 405 estimates the 

10 instantaneous phase noise by removing the linear instantaneous 
phase from the instantaneous phase. Finally, in Step 506, the 
zero-crossing resampler 406 re-samples the instantaneous phase 
noise data that is close to the zero-crossing timings of the real 
part of the analytic signal so as to estimate the timing jitter 

15 sequence, thereby finishing the procedure. 

In Step 504 where the ideal edge timing of the input clock 
signal is obtained, the initial, phase estimator 404 obtains the 
ideal edge timing in accordance with Equation (22) . In addition, 
20 Step 501 can be performed in accordance with a procedure shown 
in any of Figs. 31, 33 and 35. 



The timing estimator shown in Fig. 28 can be constructed 
as an estimator which estimates only the ideal edge timing of the 

25 input clock signal. In this case, the linear trend remover 405 
for obtaining the instantaneous phase noise by removing the linear 
instantaneous phase from the instantaneous phase and the 
zero-crossing resampler 4 06 for obtaining the timing jitter 
sequence of the input signal from the instantaneous phase noise 

30 can be omitted. Similarly, the timing estimating method shown 
in Fig. 29 may estimate only the ideal edge timing of the input 
clock signal. In this case, Step 505 for estimating the 
instantaneous phase noise by removing the linear instantaneous 
phase from the instantaneous phase and Step 506 for estimating 
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the timing jitter sequence from the instantaneous phase noise can 
be omitted. 



Fig. 3 0 shows an exemplary structure of the analytic 
5 transformer included in the timing estimator 4 00 of the present 
embodiment. An analytic signal transformer 600 includes a 
band-pass filter 601 for extracting, from the signal to be measured, 
only the components around the fundamental frequency of the signal 
to be measured, thereby band-limiting the signal to be measured; 

10 andaHilbert transformer 602 for performing Hi lbert transformation 
for the output signal of the band-pass filter 601 so as to generate 
Hilbert transform pairs of the input signal . The band-pass filter 
601 may be an analog filter or a digital filter, and may be 
implemented by using digital signal processing such as FFT. In 

15 addition, the band-pass filter 601 may be constructed in such a 
manner that the passband of the signal can be changed freely. 

Next, an operation of the analytic signal transformer 600 
of the present embodiment for transforming the signal to be measured 
20 into the analytic signal that has been band-limited is described. 
Fig. 31 shows an exemplary procedure of the signal transformation 
of the present embodiment . First, in Step 701, the band-pass filter 

601 extracts from the signal to be measured the frequency components 
around the fundamental frequency thereof, thereby limiting the 

25 band of the signal. Then, in Step 702, the Hilbert transformer 

602 applies Hilbert transformation to the signal to be measured, 
that has been subjected to band limiting, thereby generating the 
Hilbert transform pairs of the input signal, that correspond to 
the imaginary part of the analytic signal. Finally, in Step 703, 

30 the analytic signal transformer 600 outputs the output signal of 
the band-pass filter 601 as the real part of the analytic signal 
and the output signal of the Hilbert transformer 702 as the imaginary 
part of the analytic signal, thereby finishing the procedure. 
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Fig. 32 shows another exemplary structure of the analytic 
signal transformer included in the timing estimator 400 of the 
present embodiment . The analytic signal transformer 800 includes 
a time -domain to frequency-domain transformer 801 for transforming 
the signal to be measured into two-sided spectra in the frequency 
domain; a bandwidth limiter 8 02 for extracting only frequency 
components around the positive fundamental frequency of the 
two-sided spectra in the frequency domain; and a frequency- domain 
to time-domain transformer 8 03 for performing inverse 
transformation for the output of the band-litter 802 into a signal 
in the time domain. The time -domain to frequency- domain 
transformer 801 and the frequency-domain to time-domain 
transformer 803 may be implemented by using FFT and inverse FFT, 
respectively. Moreover, the bandwidth limiter 802 maybe arranged 
in such a manner that the passband of the signal can be changed 
freely. 

Next, an operation of the analytic signal transformer 800 
of the present embodiment for transforming the signal to be measured 
20 into the analytic signal that has been band- limited is described. 
Fig. 33 shows another exemplary procedure of the signal 
transformation of the present embodiment. First, in Step 901, 
the time-domain to frequency- domain transformer 801 performs FFT 
for the signal to be measured, thereby transforming the signal 
25 in the time domain into the two-sided spectra in the frequency 
domain. Then, the bandwidth limiter 8 02 replaces the negative 
frequency components of the two-sided spectra in the frequency 
domain with zeros in Step 902, and then replaces, in the one-sided 
spectra in which the negative frequency components have been 
30 replaced with zeros, the frequency components other than the 
frequency components around the fundamental frequency of the signal 
to be measured with zeros so as to leave only the components around 
the fundamental frequency in Step 903, thereby the signal in the 
frequency domain is band-limited. Finally, in Step 904, the 
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frequency- domain to time-domain transformer 803 performs inverse 
FFT for the one-sided spectra that has been band-limited so that 
the frequency-domain signal is transformed into the time-domain 
signal, thereby finishing the procedure . In this procedure, the 
5 order in which Steps 9 02 and 903 are performed may be changed so 
that Step 902 follows Step 903. Moreover, in Steps 902 and 903, 
the negative frequency components in the two-sided spectra may 
be replaced with zeros after the signal to be measured is 
band-limited by replacing the frequency components other than the 
10 components around the fundamental frequency of the signal to be 
measured with zeros so as to leave only the frequency components 
around the fundamental frequency. 

Fig. 34 shows another exemplary structure of the analytic 

15 signal transformer included in the timing estimator 4 00 of the 
present embodiment. The analytic signal transformer 1000 
includes a buffer memory 1001 for storing the signal to be measured; 
a waveform data selector 1002 for successively extracting a section 
of the signal from the buffer memory 1001 in such a manner that 

20 the section of the signal extracted presently partially overlaps 
the previously extracted section; a window function multiplier 
1003 for multiplying the extracted signal section by a window 
function; a time-domain to frequency-domain transformer 1004 for 
transforming the signal section after multiplication by the window 

25 function into the two-sided spectra in the frequency domain; a 
band-width limiter 1005 for extracting the frequency components 
around the positive fundamental frequency of the signal to be 
measured from the two-sided spectra; a frequency-domain to 
time-domain transformer 1006 for inversely transforming the output 

30 of the bandwidth limiter 1005 to a time-domain signal; and an 
inverse -window function multiplier 1007 for obtaining the analytic 
signal that has been band-limited by multiplying the time-domain 
signal by the reciprocal of the window function. The time-domain 
to frequency-domain .transformer 1004 and the frequency-domain to 
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time-domain transformer 1006 may be implemented by using FFT and 
inverse FFT, respectively. Moreover, the bandwidth limiter 1005 
may be arranged in such a manner that the passband of the signal 
can be changed freely. 



Next, an operation of the analytic signal transformer 1000 
of the present embodiment for transforming the signal to be measured 
into the analytic signal that has been band- limited is described. 
Fig. 35 shows another exemplary procedure of the signal 

10 transformation of the present embodiment. First, in Step 1101, 
the signal to be measured is stored in the buffer memory 1001. 
Then, inStepll02, the waveform data selector 1002 selects a section 
of the signal stored in the buffer memory 1001 and extracts it. 
In Step 1103, the window function multiplier 1003 multiplies the 

15 extracted section by the window function. In Step 1104, the 
time-domain to frequency -domain transformer 1004 performs FFT for 
the section of the signal, i.e., the partial signal, after being 
multiplied by the window function, thereby transforming the 
time-domain signal into the two-sided spectra in the frequency 

20 domain. Then, the bandwidth limiter 1005 replaces the negative 
frequency components of the two-sided spectra in the frequency 
domain with zero in Step 1105 and thereafter replaces the frequency 
components other than those around the fundamental frequency of 
the signal to be measured with zeros so as to leave only the 

25 components around the fundamental frequency, thereby 

band-limiting the frequency-domain signal. Then, in Step 1107, 
the frequency- domain to time-domain transformer 1006 performs 
inverse FFT for the one-sided spectra in the frequency domain after 
being band-limited, so as to transform the frequency- domain signal 

30 into the time-domain signal. In Step 1108, the inverse window 
function multiplier 1007 multiplies the time-domain signal after 
being subjected to the inverse transformation, by the inverse of 
the window function used in the multiplication in Step 1103. As 
a result , the analytic signal that has been band- limited is obtained . 
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Finally, in Step 1109, it is checked whether or not the data that 
has not been processed is stored in the buffer memory 1001. If 
the data that has not been processed is stored, the waveform data 
selector 1002 selects and extracts a section of the data 
5 successively in such a manner that the presently selected section 
partially overlaps the previously selected section, and thereafter 
Steps 1103, 1104, 1105, 1106, 1107, 1108 and 1109 are repeated. 
In a case of no data that has not been processed, the procedure 
is finished. In this procedure, the order in which Steps 1105 

10 and 1106 are performed may be changed so that Step 1105 follows 
Step 110 6. Moreover, in Steps 1105 and 1106 , the negative frequency 
components in the two-sided spectra may be replaced with zeros 
after the signal to be measured is band- limited by replacing the 
frequency components other than the components around the 

15 fundamental frequency of the signal to be measured with zeros so 
as to leave only the frequency components around the fundamental 
frequency. 

Fig. 3 6 shows another exemplary structure of the clock skew 
20 measuring apparatus of the present embodiment . The clock skew 
measuring apparatus 12 00 has the same structure as that shown in 
Fig. 25 except that analog-to-digital (AD) converters 1201a and 
12 01b are provided for digitizing an analog signal to be measured 
so as to convert it to a digital signal. (For simplification, 
25 the description for the same components as those shown in Fig. 
25 is omitted.) As the AD converter, it is desirable to use a 
high-speed AD converter, a digitizer, and a digital-sampling 
oscilloscope . 



30 Next, a clock skew measuring operation of the clock skew 

measuring apparatus 12 0 0 of the present embodiment is described. 
Fig. 37 shows another procedure of the clock skew measuring method 
of the present embodiment. This clock skew measuring method is 
the same as that shown in Fig. 26 except that Steps 1301 and 1302 
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are provided in which the analog reference signal and the analog 
signal to be measured are subjected to sampling (i.e. , digitized) 
to be converted into the digital reference signal and the digital 
signal to be measured by the AD converters 1201a and 1201b. (For 
5 simplification, the description for the same steps as those shown 
in Fig. 26 is omitted.) 

The steps for respectively converting the analog signals 
to the digital signals may be included in the procedure of the 
10 clock skew measuring method shown in Fig. 27, as shown in Fig. 



Fig. 3 9 shows still another exemplary structure of the clock 
skew measuring apparatus of the present embodiment . The clock 
15 skew measuring apparatus 14 00 has the same structure as that shown 
in Fig . 2 5 except that waveform clippers 14 01a and 14 01b are provided 
for removing the AM components of the signal . (For simplification, 
the description for the same components as those shown in Fig. 
25 is omitted. ) 



Next, a clock skew measuring operation of the clock skew 
measuring apparatus 1400 of the present embodiment is described. 
Fig. 4 0 shows another procedure of the clock skew measuring method 
of the present embodiment. This clock skew measuring method is 
25 the same as the jitter measuring method shown in Fig. 2 6 except 
that Steps 1501 and 1502 are provided in which the waveform clippers 
14 01a and 14 01b remove the AM components of the signal to be measured . 
(For simplification, the description for the same steps as those 
shown in Fig. 26 is omitted.) 



The steps for removing the signal to be measured may be 
included in the procedure of the clock skew measuring method shown 
in Fig. 27, as shown in Fig. 41. 



38 . 



20 



30 



AD-0318 



JL MJ « .3 .^^yC! 



~fl- iu» -O .rn . J2» JL<p ,io « ilji ^ji as: u c 



Fig. 42 shows another exemplary structure of the timing 
estimator of the present embodiment. The timing estimator 1600 
has the same structure as that shown in Fig. 2 8 except that a 
low-frequency component remover 1601 is provided for inputting 
5 the instantaneous phase noise and removing low frequency components 
of the instantaneous phase noise so as to output the instantaneous 
phase noise after low frequency components have been removed to 
the zero-crossing resampler. (For simplification, the 
description for the same components as those shown in Fig. 28 is 
10 omitted.) 



Next, an operation of the timing est imator 1600 of the present 
embodiment for estimating the edge timing of the clock signal to 
be measured is described. Fig. 4 3 shows another procedure of the 

15 timing estimation of the present embodiment. This timing 

estimation is the same as that shown in Fig. 2 9 except that Step 
1701 is provided for estimating the instantaneous phase noise and 
then removing the low frequency components of the instantaneous 
phase noise by means of the low-frequency component remover 1601. 

20 (For simplification, the description for the same steps as those 
shown in Fig. 2 9 is omitted.) 

Fig. 44 shows an exemplary arrangement of a clock skew testing 
system using the clock skew measuring apparatus of the present 

25 embodiment . The clock skew testing system 1800 includes automatic 
test equipment 1802 for supplying signals to a device under test 
(DUT) 1801, and the clock skew measuring apparatus 1803 for 
measuring the clock skews between the clock signals output from 
the device under test . The device under test has the clock signal 

30 selector for selecting and outputting the on-chip distributed 
clocks. The automatic test equipment • 1802 supplies a 
low-f requency systemclock f or drivinga clockdistributioncircuit , 
a control signal for controlling operation modes of the device 
under test 1801, and a test signal to the device under test 1801, 
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and also supplies the reference signal to the clock skew measuring 
apparatus 1803. It is desirable that the system clock and the 
reference signal respectively supplied to the device under test 
1801 and the clock skew measuring apparatus 1803 contain little 
jitter. The clock skew measuring apparatus 18 03 may be 
incorporated into the automatic test equipment 1802 or be 
constructed by an oscilloscope or other measuring instrument. 

As is apparent from the above, according to the clock skew 
measuring apparatus and method of the present invention, the number 
of high-frequency clock output pins that may increase the cost 
can be made minimum by providing the clock signal selector on the 
chip, thereby the cost of the clock skew test can be greatly reduced . 

Although the present invention has been described by way 
of exemplary embodiments , it should be understood that many changes 
and substitutions may be made by those skilled in the art without 
departing from the spirit and the scope of the present invention 
which is defined only by the appended claims. 



AD-0318 



